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We study theoretically the squeezing spectrum and second-order correlation function of the output light for an 
optomechanical system in which a mechanical oscillator modulates the cavity linewidth (dissipative coupling). 

We find strong squeezing coinciding with the normal-mode frequencies of the linearized system. In contrast to 
dispersive coupling, squeezing is possible in the resolved-sideband limit simultaneously with sideband cooling. 

The second-order correlation function shows damped oscillations, whose properties are given by the mechanical¬ 
like, the optical-like normal mode, or both, and can be below shot-noise level at finite times. 


I. INTRODUCTION 


II. MODEL 


Cavity optomechanics is an exciting, novel experimental 
platform that will allow us to explore fundamental questions 
of quantum mechanics and paves the way for applications in 
quantum-information processing, high-precision metrology, 
and gravitational-wave detection [ 1 ]. 

It was recognized early on that cavity optomechanical sys¬ 
tems can squeeze light similar to a nonlinear Kerr medium 
[2, 3]. Experimentally, this ponderomotive squeezing has re¬ 
cently been demonstrated as well [4-6] . As far as correlations 
between photons is concerned, photon antibunching has been 
predicted [7, 8 ] to occur in the single-photon strong-coupling 
regime [7, 9] and two-mode optomechanical systems [10-13]. 
Moreover, two-time photon correlation functions have been 
proposed as a means to observe the onset of parametric insta¬ 
bility [14] and unconventional photon blockade [12]. 

So far, most research in the field of optomechanics has fo¬ 
cused on dispersive coupling, where mechanical motion mod¬ 
ulates the resonance frequency of the cavity. Elste et al. [15] 
have proposed a novel kind of optomechanics, where mechan¬ 
ical motion modulates the linewidth of the cavity. In this case, 
the radiation pressure force spectrum features a Eano reso¬ 
nance modifying the interaction between light and mechanics 
dramatically. Cooling of the mechanical oscillator with dissi¬ 
pative coupling has been predicted [15-18] and recently also 
been demonstrated [19]. However, to date most of the proper¬ 
ties of dissipative optomechanics remain unexplored. In par¬ 
ticular, the experimental progress in this direction motivates 
us to investigate the photonic properties of this setup. 

In this work, we study theoretically the potential of dissipa¬ 
tive optomechanics as a source for squeezed light and corre¬ 
lated photons. Our analysis reveals that the system can gen¬ 
erate strong squeezing of light and we predict oscillations of 
the photon correlation function with a suppression below shot 
noise at finite times. It turns out that both the squeezing spec¬ 
trum and the photon correlation function can be understood in 
terms of the normal modes of the linearized system. 

The remainder of this paper is organized as follows. In 
Sec. II we introduce the Hamiltonian of an optomechanical 
system with dispersive as well as dissipative coupling and de¬ 
rive the linearized equations of motion describing this system. 
In Sec. Ill we present our results on the squeezing spectrum 
and in Sec. IV we discuss the properties of the second-order 
photon correlation function. We conclude in Sec. V. 


We study an optomechanical system in which a mechanical 
degree of freedom modulates the resonance frequency ooc (dis¬ 
persive coupling) and the linewidth k, (dissipative coupling) of 
a cavity mode. The Hamiltonian (h = 1) is given by [15] 

H = uJqQj^Q j ujrfjib^b H “b 
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where the dispersive and dissipative coupling strengths are 
A = —^ and B = ^, respectively. The first 

and second term in (1) describe the cavity mode (frequency 
ujc) with bosonic field operator d and the mechanical oscil¬ 
lator (frequency ujm) with bosonic field operator b. The me¬ 
chanical displacement x = xo{b ^ b^) has zero-point fluctua¬ 
tions xq = with mass m. and describe 

the damping of the cavity due to the optical and mechanical 
baths, respectively. Here, bq are bosonic fleld operators of 
the optical bath, p is the density of states of the optical bath, 
n and 7 are optical and mechanical damping rates, respec¬ 
tively. Eluctuations in the input from the optical (mechan¬ 
ical) bath are described by operators ^ and 77, respectively. 
In our analysis, we assume Markovian baths, where the me¬ 
chanical bath is characterized by a thermal phonon number 
n = with inverse temperature [3, whereas the 

optical bath is assumed to be at zero temperature. The non¬ 
zero expectation values hence are {fj^{t)fj{t')) = n S{t — t') 
and {t)r]{t')) = (n + 1 ) S{t — t') for the mechanical input, 
and {i{t) = S{t — t') for the optical input. 

Writing d = (d + with d real, we decompose 

the optical fleld into mean amplitude d and fluctuations d and 
move to a frame rotating at the laser frequency ujl- Then, 
we employ standard input-output theory [ 20 ] and linearize the 
Langevin equations of motion [15] 

^=Mu + Fui„ (2) 

dt 

where 


(d\ 


(A 


, Uin = 


c 


f) 

VcV 

Uv 











2 


and 


M = 



/ Vk 

0 

0 


F — 

0 

y/d 

0 

0 

J: — — 

y/BBa/2 

— Ba/2 

R 

0 


\—^/B Ba/2 

^/RBd/2 

0 

77/ 

/iA — k,/2 0 

El 


E 


0 -iA-Kl2 El El 

-E 2 E 2 -iuJm - 7/2 0 

V -E'2 —E 2 0 i<^m ~ 


(4) 


\ 

7 / 2 / 

(5) 


with the coefficients Ei = iAna — ^ (iA k,/2) and E 2 = 
iAK,a — ^ {iA — /^/2) and A = ujl— ojc the laser detuning. 
The linearized input-output relation for the optical field is [15] 
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( 6 ) 


From the exact solution [17] to the linearized equations of 
motion (2) we find that the following two steady-state correla¬ 
tion functions S{uj) = (dour( t) dour (0)} as well 

as N(u;) = dre^^^ (dl,^rj.(T)douT(0)} are given by 
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with the cavity response function xc'(u;) = [|—i(ci;-^A)] 
the mechanical response the op¬ 

tomechanical self-energy S(u;) = S^(u;) + S^(co’) + S^^(ci;) 
where T.a{^) = -i{AK.\d\Y[xc{<^) - Xc(-u;)], = 

i{B\a\/2f[Xc{^){iA + n/2f - Xc(-^)(^A - n/2f] , 
= BAB\a\^[xc{^){i^ + /^/ 2 ) - Xc(“^)(^^ “ 
f^/2)], and a{uj) = aA{oo) + o:b{^) with aA{oo) = 
iAf^xd^), OLBd) = f [1 -Xc(^)(^AH-/^/2 )] andQ(cj) = 
Xm(^)“^Xm(-^)“^ + with QdY = Q(-u;). 


III. SQUEEZING SPECTRUM 

In this section we will investigate the squeezing of the out¬ 
going light by a dissipative optomechanical system (OMS). 
Squeezing is conveniently characterized by the spectrum 
of squeezing Sed) = ^{AXe{uj)AXe{uj')) where 

(AX,(c^)AX,(c^')) = (A,(c^)X,(cc')) - (X,(cc))(X,(c^')) 

and Xq{uj) = + dQ^r^{uj)e~'^^^‘^. In the follow¬ 

ing we concentrate on the optimal squeezing spectrum [20] 

5'opt(<^) = minS'0(co’) = 1 — 2 |5'(co’)|+A(co’)+A(—cc) (8) 
9 

which can be expressed in terms of the two correlation func¬ 
tions S{uj) and N{uj) given in Eq. (7) above. 


Choosing the mechanical frequency ujrn to be unity, the op¬ 
timal squeezing spectrum 5 'opt(<^) is a function of detuning 
A, cavity linewidth b, coupling strength Aa and Ba, as well 
as thermal phonon number n. Note that the single-photon cou¬ 
pling strengths, A and B, and the intracavity amplitude a only 
appear as a product and not individually in Eqs. (7) and (8). 

Fig. 1 (a) shows the optimal squeezing spectrum 5 'opt(<^) 
for a dissipative OMS as a function of detuning A. Shaded ar¬ 
eas correspond to regions where the solution of the linearized 
equations of motion (2) is unstable, i.e. the eigenvalues of the 
matrix M (5) do not all have negative real parts. Focusing on 
the stable regions, we find two different types of behavior. On 
the red-detuned side A < 0 we observe two dips below shot- 
noise level 5 'opt(<^) = 1. In Fig. 1 we also plot the normal¬ 
mode frequencies of the system obtained from the eigenvalues 
of the matrix M (5). We see that the squeezing dips coincide 
with the normal-mode frequencies of the system. This is qual¬ 
itatively similar for the dispersive OMS shown in Fig. 1 (b). 
Mathematically, this is a consequence of the fact that Q(u;) is 
the determinant of the matrix M. Physically, this corresponds 
to the fact that the OMS will respond most strongly close to 
its resonances, i.e. normal-mode frequencies, and the large re¬ 
sponse can in turn strongly affect the outgoing light field. 

In Fig. 1 (b) we see that in the limit uOm > k, squeezing is 
small for dispersive coupling at A —^m- Stronger squeez¬ 
ing can be achieved around resonance A ^ 0 close to the 
mechanical frequency 00 ^ ujm and in particular close to zero 
frequency u; « 0 [2, 3]. For dissipative coupling, and in con¬ 
trast to the dispersive case, there is a window of stability on 
the blue-detuned side A > 0. For small positive detuning A 
strong squeezing occurs at frequencies close to the mechani¬ 
cal frequency uj ujm- A special point is A = ci;^/2 where 
the squeezing spectrum has a single dip close to the frequency 
of the optical-like mode cc « A. Moreover, strong squeezing 
also occurs at the point where the two normal-mode frequen¬ 
cies merge close to the onset of the parametric instability. 

Fig. 1 (c) shows the optimal squeezing spectrum 5 'opt(<^) as 
a function of the dissipative coupling strength 5a for a system 
driven on the red sideband A = —At weak coupling the 
two normal modes are degenerate leading to a single dip in 
the squeezing spectrum. At larger coupling the degeneracy 
is broken and as normal-mode splitting (NMS) develops two 
dips emerge in the squeezing spectrum 5 'opt(<^)- NMS in the 
squeezing spectrum is also present in the dispersive case, see 
Fig. 1 (d). Note that squeezing on the sideband A = —uOm is 
much weaker in general for dispersive coupling. 

Fig. 1 (e) shows the optimal squeezing spectrum 5 'opt(<^) 
of a dissipative OMS as a function of cavity linewidth b. For 
B ujjn we find two narrow dips in the squeezing spectrum. 
For larger cavity linewidth n the dips broaden and merge into a 
single broad region of squeezing. This is in contrast to the case 
of a dispersive OMS, see Fig. 1 (f), where there is only weak 
squeezing in the resolved sideband limit n ujm, and a broad 
dip in the squeezing spectrum occurs in the bad-cavity limit 
B ^ uJrn- We note that since ground-state cooling for disper¬ 
sive coupling is possible only in the resolved-sideband limit, 
a dissipative OMS offers the advantage that strong squeezing 
and ground-state cooling can be realized simultaneously. 
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FIG. 1. Optimal squeezing spectrum 5 'opt(^) as a function of detuning A/ujm for n/ujm = 0.3, Ba = 0.6 (a) and Aa = 0.6 (b), as a function 
of coupling strength Ba (c) and Aa (d) for K,/oJm = 0.3 and Ajujrn = — 1, and as a function of cavity linewidth n/um for A/c^m = — 1 and 
Ba = 0.4 (e) and Aa — 0.4 (f). Throughout this work, we consider a mechanical resonator with high quality factor — 10“^ and at 

zero temperature n = 0. The shaded areas correspond to regions of instability. The black lines indicate the normal-mode frequencies. 


IV. PHOTON CORRELATIONS 

We study next the photon correlations in the outgoing light 
of a dissipative OMS. Applying Wick’s theorem, which is ex¬ 
act for the Gaussian states considered here, and using the de¬ 
composition of field operators a = + one obtains 

the second-order photon correlation function g^‘^\r) [12, 21] 

= 1 I 2|«0UT|"Re [N{t)] + 2Re [(a*)SuT5(^)] 

^ [|aouT|2 + 7V(0)]2 

, |A^(r)P + |5(r)|^ 

[|aouT|2 + 7V(0)]2 

where the output field given by aoux = f + |), 5 '(t) = 

(<^out(^)<^out( 0)), and N{r) = {d^ujir )(iouT(0))- The cor¬ 
relation functions can be found by Fourier transforming the 
analytic expressions (7) or solving numerically their equation 
of motion that can be derived from Eq. (2). 

We find that g^‘^\r) shows damped oscillations with one 
or several frequencies depending on the choice of parameters, 
see Fig. 2 (a), (c), (e). To identify the spectral contributions 
we calculate numerically g^‘^\uj), the Fourier transform of 


g^‘^^ (r) — 1, as shown in Fig. 2 (b), (d), (f). We see that in (a) 
the dynamics of g^‘^\r) is dominated by one spectral com¬ 
ponent close to the mechanical frequency In Fig. 2 (b) 
we also plot a Lorentzian whose linewidth and frequency are 
given by the mechanical-like normal mode, and in Fig. 2 (a) 
we plot the corresponding exponential envelope. From this 
analysis, we conclude that in this case g^‘^\r) is dominated 
by the mechanical-like normal mode, i.e. its frequency is the 
(effective) mechanical frequency and the linewidth is the (ef¬ 
fective) mechanical linewidth. In contrast, in Fig. 2 (c) and 
(d), the frequency and linewidth are determined by the optical- 
like normal mode and its properties. In Fig. 2 (e) and (f), both 
normal modes contribute significantly leading to a beating in 
g^^Kr) , as it can be nicely seen in the spectral domain g^‘^\uj). 

It is worth noting that we observe g^‘^^ (r) < 1 at finite times 
r, i.e. the OMS suppresses the probability of photon arrival in 
the outgoing light beam which may be useful for applications. 

In addition to the response at the normal-mode frequencies, 
Fig. 2 (b) shows small features at zero and twice the normal¬ 
mode frequency. These are due to the non-linear terms | 5 '(t) p 
and \N{r)\‘^ in (9). While they are important at the single¬ 
photon level a ^ 1, they become gradually less important in 
the limit of large intracavity field a 1 [21]. 
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FIG. 2. Photon correlation function as a function of time r (blue solid) for K./uOm — 0.8, A/um = 0.1, Ba = 0.4 (a), hv/um = 0.3, 

Ajujrn = 0.5, Ba = 0.9 (c) and K./uJm = 0.3, Ajurn = 0.3, Ba — 0.9 (e). Orange dashed lines are Ce~^^ where F is the deeay rate of the 
mechanical-like (a) or optical-like (c) normal mode. Panels (b, d, f) show the correlation spectra (cc) (blue solid) eorresponding to g^‘^^ (r) 
in (a, c, e). Red dashed lines show a Lorentzian with the frequency and linewidth of the mechanical-like (b) or optical-like (d) normal mode. 
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FIG. 3. Photon correlation funetion g^‘^\r) (blue solid) and approximation g^‘^\r) (orange dashed) as a function of time r for (a) d = 1.0, 
(e) a = 2.5, and (e) a = 10.0. The correlation spectrum g^‘^\uj) (blue solid) and approximation g^‘^\(jj) (red dashed) in (b), (d), (e). The 
parameters are A^/cCm = 1-2, A/um = —1, Ba = 0.4. While non-linear effects are strong for d = 1 (a), they are negligible for d = 10 (e). 


To investigate their effect further we plot in Fig. 3 (a), (c), 
(e) the exact photon correlation function g^^\r) as well as 

-Bl, > . , 2|aouTpRe[W(T)| + 2Re[(a*)?,„TS(T)] 

+-PrfTw(0)P- 

( 10 ) 

neglecting the nonlinear terms |5 '(t)P and |W(r)p in (9). In 
Fig. 3 (b), (d), (f) we plot the exact correlation spectra g^‘^\uj) 


as well as g^‘^^ (uj), the Fourier transform of g^‘^^ (r) — 1, 

" [|aouTP+iV(0)]^ ■ 

( 11 ) 

In Fig. 3 (b) we clearly observe the non-linear contributions, 
|5'(r)p and |W(r)p, at zero and twice the normal-mode fre¬ 
quency which modify g^^^ (r) dramatically, see Fig. 3 (a). 
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FIG. 4. Photon correlation spectrum for dissipative OMS as a function of (a) detuning A for K./uJm — 0.1, Ba = 0.4; (c) as a 

function of coupling strength Ba for K./uJm = 0.3, A/ujm = — 1, and (e) as a function of cavity line width k for A/cum = —1, Ba = 0.4. 
Corresponding correlation spectrum (cc) for dispersive OMS in (b), (d), (f). We consider a high quality factor 7/ct;m = 10“^ and normalize 
g^‘^^ (cc) to unity in each of the panels. Shaded areas correspond to regions of instability, and black lines indicate the normal-mode frequencies. 


We note in passing that for single-photon coupling strength 
B = 0.4 as shown in (a) one should compare our results to 
numerical solutions of the quantum master equation. 

In contrast to the squeezing spectrum 5 'opt(<^) (8 ), the pho¬ 
ton correlation function g^‘^^ (r) depends not only on the prod¬ 
uct Ba of coupling strength and intracavity amplitude, but 
also on the intracavity amplitude a itself. Inspecting the scal¬ 
ing in (9) we find that photon correlations can be large for 
small cavity amplitudes a ^ 1 , but approach those of coherent 
states g^‘^^ (r) —> 1 for a ^ oo with Ba fixed, see Fig. 3 [21]. 
For a large intracavity field a 1, where g^‘^\uj) is a good 
approximation, the photon correlation function g^‘^\r) only 
depends on detuning A, cavity linewidth k, coupling strength 
Ba, and thermal phonon number n. 

In Fig. 4 (a) we show the correlation spectrum g^‘^^ (cc) of a 
dissipative OMS as a function of detuning A. The shaded ar¬ 
eas correspond to regions where the solution to the linearised 
equations of motion (2) is unstable. Similar to the squeezing 
spectrum 5 'opt(<^) the correlation spectrum g^‘^^ (cc) can be un¬ 
derstood in terms of the normal modes of the system whose 
frequencies we also plot in Fig. 4. We identify regions where 
the (strongly broadened) mechanical-like normal mode domi¬ 
nates the dynamics of the photon correlation dynamics around 


A ^ 0, like in Fig. 2 (a), the special point A = cc ;^/2 where 
the optical-like normal mode dominates, like in Fig. 2 (b), or 
other cases, like in Fig. 2 (c), where both modes contribute 
significantly. 

Before concluding we note that as expected with increasing 
thermal phonon number n the quantum features in the second- 
order correlation function as well as the squeezing spectrum 
gradually vanish. 


V. CONCLUSION 

We investigated dissipative OMS as a source of squeezed 
light and correlated photons. Strong squeezing occurs at the 
normal-mode frequencies of the system, and in contrast to dis¬ 
persive OMS, even in the sideband-resolved limit. The pho¬ 
ton correlation function g^‘^\r) exhibits damped oscillations 
whose frequency and linewidth are given by the properties of 
the normal modes. Depending on parameters this is the me¬ 
chanical frequency and its (effective) decay time, or the op¬ 
tical detuning and cavity linewidth, or both modes contribute 
significantly leading to a beating in the correlation function. 
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